Fano Resonance Between Mie and Bragg Scattering in Photonic Crystals 
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We report the observation of a Fano resonance between continuum Mie scattering and a narrow Bragg band 
in synthetic opal photonic crystals. The resonance leads to a transmission spectrum exhibiting a Bragg dip with 
an asymmetric profile, which can be tunably reversed to a Bragg rise. The Fano asymmetry parameter is linked 
with the dielectric contrast between the permittivity of the filler and the specific value determined by the opal 
matrix. The existence of the Fano resonance is directly related to disorder due to non-uniformity of a-Si02 opal 
spheres. Proposed theoretical "quasi-3D" model produces results in excellent agreement with the experimental 
data. 

PACS numbers: 42.70Qs, 42.25.Fx, 42.79.Fm 



Mie and Bragg scattering are key optical phenomena in 
photonic crystals (PhC) composed of spherical or nearly 
spherical particles. Light scattering by an isolated spherical 
particle can be described by Mie theory [ 1 ] . Considering PhC 
composed of a periodic array of such spheres, interference of 
scattered waves results in the transformation of Mie scatter- 
ing into Bragg scattering and gives rise to formation of the 
photonic band structure [2]. The underlying Mie scattering is 
therefore hidden in perfectly ordered PhC and as a result, it 
has been insufficiently studied so that its role is clear only for 
the case of perfect structures. The resulting Bragg scattering, 
on the other hand, has been intensively examined, both exper- 
imentally and theoretically in great detail [IH El H 0] ■ In 
particular, the multiple Bragg diffraction phenomenon in PhC, 
which occurs when two narrow Bragg bands demonstrate the 
avoided crossing effect [7,^^ JDJ has been thoroughly stud- 
ied. Departing from this phenomenon and with intention to 
further deepen our understanding of light scattering in PhC, 
a number of challenging problems can be formulated: What 
can we expect if a spectrally naiTow Bragg band interacts with 
a broad spectrum originating from certain scattering mech- 
anisms such as Mie or Fabry-Perot scattering? Is it possi- 
ble to observe the consequences of this interaction or simply 
Mie scattering experimentally? What are the effects of inher- 
ent disorder in the structural components of opal-based PhC, 
beyond the well-known broadening and degradation of stop 
bands fl El El El El El ? 

If a narrow Bragg band interacts with the continuum spec- 
trum through an interference effect constructively or destruc- 
tively, we can expect an interaction of Fano-type [16], a 
phenomenon well-known across many different branches of 
physics. The Fano resonance between continuum and discrete 
states manifests as an asymmetric profile of narrow band in 
the transmission spectrum, which in general has the form: 



(1) 



where Vl — {ix) — oj b) / {jb /2), ojb the frequency, jb is the 



width of the narrow band, and q is the Fano asymmetry param- 
eter. It was shown theoretically that Fano-type asymmetric 
line shapes can be created in the response function of certain 
PhC. In general these systems consist of a waveguide with for- 
ward and backward propagating waves being indirectly cou- 
pled via one or more mediating resonant cavities, defects or 
another similar mechanism | I4EI. S 23. 21, 22, 23]. Sur- 
prisingly, despite this growing interest in theoretical and ex- 
perimental studies of Fano resonances in PhC, until now no 
one has ever considered a Bragg band acting as the discrete 
narrow band, as opposed to resonant cavities or defect states. 

In this Letter, we report on the experimental discovery of 
Fano resonances involving interference between Mie scatter- 
ing (as a continuum) and Bragg scattering (as a discrete 
state). The (111) Bragg band in the transmission spectra 
of synthetic opals demonstrates all characteristic features of 
the Fano resonance, including asymmetry of the transmission 
band which in this case can be inverted around its central fre- 
quency when the sign of the Fano parameter q is reversed. 
Moreover, we report observation of an impressive effect of 
the Bragg band reversal, which manifests as appearance of a 
transmission "rise" instead of a transmission dip, and which 
we refer to as ''enhanced Bragg transmission". We also pro- 
pose a theoretical "quasi-3D" model able to reproduce all spe- 
cific features of the experimentally observed Fano resonance. 

To study the general features of light scattering in PhCs, we 
chose synthetic opals as an archetypal 3D PhC [6:^ H El EH 
El [3. 0; is^l • An fee lattice of opals is formed from nearly 
spherical a-Si02 particles which are neither uniform in size 
nor homogeneous in dielectric permittivity 1 11, 26, 27]. To 
obtain precise information about the inherent inhomogeneity 
of the fl-Si02 particles and statistical characteristics of the en- 
semble of particles in the investigated samples, scanning elec- 
tron microscopy (SEM) micrographs were analyzed using an 
original procedure based on the Hough transform |28]. From 
the micrograph shown in Fig. 1(b) we obtain several impor- 
tant structural parameters of the sample under study, such as 
the average diameter of the a-Si02 particles Z3=316.2nm and 
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FIG. 1: (Color online) (a) The SEM image of the (111) opal plane 
with a-Si02 particles, (b) The circles and the centers of circles de- 
rived by processing the SEM image (a) using Hough-type procedure, 
(c) The distributions of the circle diameters D (red circles) and the 
distances between circle centres ooo (blue squares) obtained by pro- 
cessing the image (b). The curves are the results of fitting with a 
Gaussian function. 



full-width at half-maximum of its normal (Gaussian) distri- 
bution of 22.8 nm. The latter corresponds to deviation in the 
diameter D equal to 7.2% at half maximum. The average dis- 
tance between particle centres aoo=315.2 nm, while D > aoo 
because of sintering of neighbouring fl-Si02 spheres in opals. 

We investigated transmission spectra of two samples of a- 
Si02 particles with average diameters of 260 nm and 316 
nm respectively. The sample was placed at the center of a 
vessel with plane-parallel quartz windows, which was filled 
with an immersion liquid. The liquid served as the opal filler 
and as the ambient medium for the sample. The permittivi- 
ties of the liquids were precisely measured by an Abbe IRF- 
454B2M refractometer 1 1 1, 26]. Two liquids - distilled water 
with e /=1 .778 and propylene glycol with e /=2.053, as well as 
their mixtures were used, providing variation in the dielectric 
constant of the filler in the range of 1.778 < £/ < 2.053. 

The transmission spectra of synthetic opals were recorded 
by a Perkin-Elmer Lambda-650 spectrophotometer in the 
215-850 nm wavelength range when white light was passed 
through an oriented opal sample along the [lll]-axis. In 
this scattering geometry, the (111) stop band is dominant 
(Fig. 2), and its spectral position is determined by ^{m) ~ 
^/8/SD./gF^+JT^^g)sJ, where g « 0.74 is the fill frac- 
tion of the fl-Si02 particles which have an average dielectric 
permittivity e~s=1.92 |11]. Careful measurements of the real 
es{r) profile function remains a challenge due to the small 
dimensions of the spheres. A possible information source 
might be the particle vibration modes I29i] . Higher Miller- 
index (hkl) stop bands are also seen in the 215-450 nm wave- 
length range [11]. Besides all these stop bands, the transmis- 
sion spectrum features a continuum broad part with an abrupt 
decrease in the transmission at about A < A(m')/2. This spe- 
cific shape was reproduced in a number of works on spec- 
troscopy of PhC composed of spherical particles and there- 
fore can be attributed as a feature typical of these structures 
Js, 13, Hi 12, 24, 25, 2^. Figure 2 shows a sharp and clear 
(111) band in the transmission spectrum, and this fact allowed 
us to analyze its shape especially thoroughly. The results 
clearly indicate a remarkable transformation of the (111) band 
when the filler permittivity is changed [Fig. 3(a)]. At first. 
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FIG. 2: (Color online) Transmission spectra of two opal samples 
(_D=260 nm, and D=316 nm) for the T ^ L scattering geometry and 
the results of calculations of the transmission spectra of an disordered 
ensemble of isolated spherical particles with _D=260 nm and D=316 
nm with use of Mie theory. 



the well-known effect of diminishing of the stop band inten- 
sity is observed when the filler permittivity approaches a cer- 
tain value, which was determined for the family of (1 1 1) stop 
bands of the 316 nm opal sample to be (111) = 1.816. The 
dependence of this value on (hkl) index of stop band family 



is due to the inhomogeneity of the fl-Si02 particles 111 



In this work, however, we focus on a different phenomenon 
which can be clearly observed in the transmission spectra, but 
has not, to our knowledge, previously been discussed. We ob- 
serve an asymmetry of the (111) dip which is inverted around 
a central frequency when the sign of the dielectric contrast 
£/ — is reversed. For (e/ — e°) < (green curves in Fig. 
3), the long-wavelength wing of the (111) stop band is rela- 
tively flat in contrast to a steep short-wavelength wing. For 
(e/ — Ey) > (blue curves in Fig. 3), the situation reverses 
- the long- wavelength wing of the (111) stop band becomes 
abrupt as compared to the now relatively flat short-wavelength 
wing of the band. Quite surprisingly, we found that there is no 
value of £/ at which the (111) band completely disappears. 
Moreover, when the filler permittivity is e/ = e° (red curves 
in Fig. 3), a transmission Bragg rise is observed instead of the 
conventional transmission Bragg dip; i.e. we have discovered 
the "enhanced Bragg transmission" . Note that the transmis- 
sion Bragg rise could only be observed in a very narrow range 
of the filler permittivity values Ae/ k, 1.816 ± 0.003. There- 
fore, this phenomenon was overlooked in our previous papers 
(Refs.[lill26D, where the filler permittivity was changed using 
a larger step (Ae/ ~ 0.01 - 0.03). 

From a theoretical point of view, light scattering in PhC 
is inherently a multi-level process. The first level represents 
scattering by a single particle with its own scattering proper- 
ties. At the higher levels, multiple scattering events by the 
identical particles arranged in a periodic structure add up, re- 
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FIG. 3: (color) (a) The transmission spectra of an opal sample 
(Z)=316 nm, thickness ~ 0.6 mm) as a function of the filler permit- 
tivity £/ in the region of the (111) photonic band (black curves) . The 
color curves (red, blue and green) are the results of fitting with a Fano 
function, (b) The results of calculations of the transmission spectra 
using the "quasi-3D" model of disordered opal structure described in 
the text: number of layers was taken to be 1000 and averaging was 
done over 40000 realizations for £/ = 1 . 8 1 6 and 1 0000 realizations for 
other values of e/ . In (a) and (b) the curves are shifted vertically by 
the values shown. 

suiting in formation of the Bloch waves of a perfect PhC. Dis- 
order, however, induces additional random phase shifts and 
amplitude variations at the first-level, single-scattering event. 
This, in turn, modifies higher level multiple scattering events 
and causes degradation of bands and Bloch modes, as well as 
giving rise to the appearance of broad background continuum. 
Finally, Bloch waves and the background continuum can in- 
terfere, giving rise to a Fano-type resonance. 

It is now clear that in order to describe the experimentally- 
observed phenomenon, a theoretical model should incorpo- 
rate both scattering mechanisms (Bragg and Mie), in a way 
that naturally permits their coupling. In the model we pro- 
pose here, the scattering properties of a disordered array of 
spherical particles forming a two-dimensional (2D) layer, cor- 
responding to the (111) plane of the fee opal structure, are 
calculated at the first stage with use of the Mie theory 
In-layer diffraction and multiple scattering processes are ne- 
glected in this model, which restricts the model applicabil- 
ity to the low contrast regime and implies that we neglect or- 
dering of the particles within the layer It can be shown that 
in this approximation, the transmittance and reflection of the 
layer are directly related to the forward and backward scatter- 
ing amplitudes S{9 = 0) and S{0 = 180°) of Mie theory 111], 
respectively, through 

tn = e'"" (l - ^S„{9 - 0)) , (2a) 



Tn = e'^''^Sn{0 = 180°). (2b) 

In the last expressions n is the layer index, k is the modulus 
of the wavevector in the filler material, h is the layer thick- 
ness, rj is the density of the particles in the layer, and the tilde 
above the scattering amplitudes implies intralayer averaging 
over the particles' size r and dielectric constant e. This aver- 
aging is conducted for every single layer by numerical inte- 
gration of the scattering amplitudes around average values e„ 
and f„ which vary from layer to layer: 

/CO /.oc 
/ S{9,e,r)f„{e-e,„r~f,,)dedr, (3) 
-oo ./ — oo 

here /„(e, r) is a probability density function taken in our cal- 
culations to be normal (Gaussian) distribution. 

At the second stage, the optical properties of the fee pho- 
tonic crystal built from a sequence of such (111) layers are 
determined by the conventional 2x2 transfer matrix technique 
ll30ll . The transfer matrix of every layer of the structure is 
expressed through transmission and reflection coefficients ob- 
tained at the previous stage with account of interlayer disor- 
der, i.e., as was explained before, the permittivity of parti- 
cles and their size also vary from layer to layer Finally, the 
calculated transmittance is averaged over a sufficient number 
of realizations. This model therefore can be referred to as 
"quasi-3D" since the Bragg diffraction is purely ID, while the 
disorder induced scattering is 3D (both intralayer and inter- 
layer disorder are taken into account). In our calculations, the 
particle size distribution was taken to be normal (Gaussian) 
with full-width at half maxima of 7.2% for average diameter 
of £'=316.2 nm (Fig.l). The same model was taken for mod- 
eling the distribution of permittivity. The only adjustable pa- 
rameter was average sphere permittivity. As can be seen from 
calculated spectra shown in Figs. (2,3), all peculiarities of the 
experimental spectra, such as asymmetry of the Bragg dip, ap- 
pearance of the rise corresponding to enhanced Bragg trans- 
mission, and the overall behavior of the transmission spectrum 
are revealed by our model. Additional calculations using other 
permittivity distributions Es (such as rectangular or triangu- 
lar), also demonstrated the presence of the Fano resonance, 
including the appearance of the Bragg rise. 

We have now elucidated a scenario of light scattering in 
realistic PhC, which are always imperfect due to the pres- 
ence of inherent disorder of their constitutive elements. The 
disorder breaks a condition of ideal Bragg scattering, which 
determines all properties of perfect structures, and gives rise 
to an additional scattering component. This defect-induced 
extra-scattering produces background radiation and is respon- 
sible for formation of a continuum spectrum with its specific 
characteristics determined by the light scattering by elements 
constituting PhC. In particular, in disordered ID PhC com- 
posed of dielectric slabs, it would be extra Fabry-Perot scatter- 
ing due to fluctuations in slabs permittivity, while in 3D PhC 
formed from spherical particles, this mechanism is Mie scat- 
tering. That is, it is the defect-induced Mie scattering which 
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is responsible for and determines the character of the contin- 
uum spectrum of opals (Fig. 2). The latter is characterized by 
the abrupt decrease in the transmission observed in the short 
wavelength part of the spectra A < Avn / 2 for various PhC 
formed from spherical particles ||l[i3,[ljltill2l|2l|2il. An 
increase in the scattering intensity responsible for this feature 
is a well-known from Mie theory fT]. This is why transmis- 
sion through collection of particles obtained within Mie the- 
ory for both £'=260nm and £'=316nm nicely reproduce this 
feature. The continuum defect-induced Mie spectrum inter- 
fere with narrow Bragg bands, and, according to our calcu- 
lations, this effect is related, in the first place, with inhomo- 
geneity of the fl-Si02 particles with respect to dielectric per- 
mittivity Es- In addition we obtained good fitting of the exper- 
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FIG. 4: (Color online) The Fano asymmetry parameter q that changes 
continuously from negative to positive values with filler permittivity 
£/ increasing for three samples with different thickness of 0.6, 0.8, 
and 2.2 mm, Z)=316nm. The solid line is a guide for the eyes only. 

imental spectra by calculations with use of Eq. (1). This fact 
confirms that Fano-type interference between the continuum 
defect-induced Mie spectrum and narrow Bragg band takes 
place. Moreover, when Fano parameter q sign is reversed 
[which corresponds to reversal in the sign of dielectric con- 
trast (Fig. 4)], this results in inversion of (111) transmission 
dip (Fig. 3). The fact that sign reversal of parameter q coin- 
cides with sign reversal of dielectric contrast {e/ — indi- 
cates that when Sf — e'j reversal of the phase difference be- 
tween disorder-induced background radiation and Bloch wave 
takes place. As it follows from analysis of curves obtained 
from Eq. (1) when q=Q narrow band should transform from 
transmission dip into transmission rise, which indeed is ob- 
served in both experimental and calculated spectra (red curves 
in Fig. 3). 

To conclude, we report a beautiful example of a Fano-type 
resonance in PhC, which originates from interference between 
narrow Bragg band and the continuum spectrum. This dis- 
tinctive continuum spectrum observed in opals is for the first 
time related with residual defect-induced Mie scattering. The 



"quasi-3D" theoretical model proposed for description of the 
light scattering in opals with account of disorder allowed us 
to depict all experimentally observed peculiarities. We found 
that the 3D nature of the disorder, especially in-plane posi- 
tional disorder and from plane to plane variation of permittiv- 
ity, is essential for explanation of the experimental results. 

In a pioneering work |31[1 on light scattering in disordered 
PCs it was predicted that localization of photons may occur in 
the vicinity of band edges. Since then such interplay between 
Bragg diffraction and disorder-induced scattering was consid- 
ered only as a way to localize light and as a mechanism ca- 
pable to suppress transmittance and increase back-scattering. 
Our studies, however, demonstrate that it is not always so 
and at some particular conditions Fano interference between 
scattered background radiation and Bloch waves gives rise to 
anomalous increase of the optical transmittance near photonic 
band edges or when Bragg condition is satisfied. 
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